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ABSTRACT 


The generally difficult problem of a system with time 
delay is studied by using parameter space methods and com- 
puter analysis. The solution is proved to be very easily 
obtained and the results presented are extendable under 
certain circumstances to other systems and useful in other 


fields of analysis. 
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1. INTRODUCTION 


Dead time or transport lag, defined as a fixed time de- 
lay in the transmission of information between two segments 
of a system, is a phenomenon that occurs frequently in every 
field of engineering. The synthesis of a control system 
where dead time tends to occur has most discouraging aspects 
and the first thoughts of the engineer confronted with the 
problem should be how to go about changing the system design 
so as to eliminate its sources. 

However, many times this is not possible due to physical 
limitations in the system. The dead time may be due to hav- 
ing the measuring element and the controller actuator neces- 
sarily separated by a certain distance, as in the extreme 
case of the remote control of a space vehicle. In other 
cases, there is a.delay between the time at which the con- 
troller actuator is moved and the time at which the plant 
reacts to the change, as in the case of an engine-governor 
system. : 

Dead time deteriorates seriously the stability of a 
system since it creates a phase lag that constantly increases 
with w. 

The analytic treatment of the problem is strenuous and 
obtaining an exact solution is very difficult, if not impos- 
Sible. The Characteristic Equation (C.E.) of a system with 


time delay has transcendental nature and a Maclaurin series 
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of the exponential term shows that the system has an infinite 
number of roots. 

Bode and Root Locus analysis are practical and efficient 
methods in the study of linear systems, but in the case of a 
system with transport lag they become too cumbersome. They 
are also single parameter techniques. However, it is pre- 
cisely in the case of systems with transport lag that the 
"Parameter Space Method" (1), (2), appears as a most useful 
and practical aid. This method has application whenever the 
system has several adjustable parameters and consists gener- 
ally in the mapping of the S plane in a N dimensional space 
of the parameters of the system. With this mapping completed 
it is possible to choose the N parameter values that allow 
the system to have a certain set of root locations. 

As previously mentioned, a system with transport lag has 
a C.E. with an infinite number of roots; therefore, if we 
were attempting to look at stability, the knowledge of the 
locations of some of the roots would not be the answer. 
Methods are available for testing the stability of systems 
with transport lag (2),(3). In this work we will concen- 
trate on finding the parameter values that guarantee a near- 
optimal behavior of the system recognizing, however, that 
analysis of stability should be performed after this synthe- 
sis phase is finished. 

The study of the general problem was followed by appli- 


cation to a practical example. 





II. THEORETICAL BACKGROUND 


Consider the two block diagrams of Figure l. The trans- 


fer functions of the two systems represented are: 


EUIS 
Т.Б.) = созше EE (1) 
1+G(S).H(S).e 
_ G(S) 
NS me (2) 


1+G(S).H(S).e 

We see that the form of the C.E. of the system does not 
Bend on the location of the dead time block, that is, the 
C.E. is the same independently of having the dead time block 
in the forward or αερα τὰ path. This allovs the extension 
of the following study to a large class of problems where 
dead time occurs. 

The term exp(-TS) gives to this equation a transcendental 
nature. In fact, the Maclaurin expansion of exp (-TS) is: 


We see then that a system with transport lag has an in- 
finite number of roots which, other than making complete sta- 
lity analysis very m it, has the effect of making some 
graphical methods, as Root Locus and Bode παρα very dus 
appropriate. 

The stability of systems where dead time appears is al- 
Ways a factor to consider since stability is seriously 


affected by it. Consider that if we examine P through 
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sinusoidal analysis and let s - jo, the transport lag term 


becomes 


e jut 


I 


F (jw) 


pu aes (4) 


We see then that the transport term has a magnitude of 
1 and a lag phase angle of WT. The phase lag added by the 
transport lag term increases as ü increases. The system 
may, therefore, be very difficult to stabilize. 

However, we assumed initially that the nature of the 
System is such that all the roots lie in the left half of 
the S plane. The system is then intrinsically stable and 
the problem is to find the values of the adjustable para- 
meters that make the system behave in an optimal manner in 
relation to an established criterion. We concentrated on 
finding the parameter values that guarantee the existence 
of a set of dominant roots for the system, that is, a set 
of roots close to the origin. Stability analysis, using for 
instance, the methods in references (2) and (3), should be 
used after the parameters for the system were selected. To 
determine the parameter values, the parameter space method 
can be used with great -advantage since, as will be shown, 
it allows us to relate precisely parameter values with root 
locations. We considered the case of systems with three 


adjustable parameters. 
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The general form of the C.E. of a system with transport lag 
is: | 

Ν Ν 
F(S) = у K,s* + у KS ce = 0 (5) 
0 j-0 


й = 
For DNE Obl.) 


οὐ τα x 
e х8 E ΝΕ mU (cos Tw + j sin Tw) (6) 


This transformation allows the separation of the C.E. 
into two equations, one containing the real and the other 
the imaginary part. That is: 

F(o,w) = F (0,0) + jF (0,0) = 0 


implies that: 


! 
© 


(7) 
(8) 


Е (О,ш) 


n 
© 


F (0,0) 


If we consider that the system has three adjustable parameters 
À, B and C, then they appear in F, (0,0), and/or F (0,0). 
Fixing one of the parameters (for instance, parameter A), 
equations (7) and (8), may be solved for the two parameters 

B and C. Recapitulating, we have that equations (7) and (8) 
are functions of the variable parameters A, B and C of O and 
W from the assumed roots -O + jw. Fixing all these variables 
but two of the parameters, for instance, B and C, simulta- 
neous solution of the two equations allows the determination 
of the two values B and C that together with the fixed 
value of A guarantees that the system has a pair of conju- 


gate roots at -O + jo . 
o — ο 
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Repeating this process for all different values of w in 
a certain range permits the determination of different set- 
pairs, B oC as in Figure 2. Of course, there is an infinite 
number of these solutions. The union of these solutions 
B.C, is a line in the parameter plane of axis B and C as 
is shown in the Figure 3 (o line). Working similarly, but 
with the variation of o instead of w, it is possible to 
obtain similar lines in the parameter plane B-C for which ἡ 
is fixed (w lines). 

The τ. of a certain line of fixed ω (ш line) 
and another line of fixed O, (o line) gives a point in the 
Parameter plane from which it is possible to read the pair 
of parameters i ONG that guarantee the existence of a pair 
of roots for the system of the form "σα Е jo,» See Figure 4. 

Repeating the above process for different values of A 
allows the determination of a line (or a group of points) 
in the three-dimensional space with ordinates A, B and C for 
which every point corresponds to assigning the system the 
pair of roots ms. + jO ° See Figure Š. 

If we assume that the system has a third dominant root 
Бо of a real nature, we can, fixing the value of A, deter- 
mine a line of fixed 6(6 | line) in the parameter plane B-C 
as in Figure 6. This is so because if we substitute the 
solution ò 2 in the characteristic equation and fix its value 
together with the value of A, then this equation becomes the 


equation of a line in the plane B-C. 


16 





Assuming different values for A, but keeping ô constant, 
the repetition of this process for each value of A results 
in a group of lines in which 6 is constant, and these lines, 
if this repetition were performed for an infinite number of 
values of A, form a surface (surface 6) in the parameter 
space A, B and C for which the system has a real root $6» 
as in Figure /7. 

In the three-dimensional parameter space, the pierce 
point(s) of the line -o, * ju, with the plane -5_ gives one 
or more points that correspond to the set(s) of parameters 


A B and Co» for which the system has a set of three roots 


o? 
po^. + jo, and б that we assumed would be the dominant ones. 

Notice that, since we are mostly interested in finding 
roots that lie close to the origin, the range of variation 
of 0, w, and 6 during the process of mapping the S-plane in 
the parameter space does not need to be very large. 

Of course, such an extensive manipulation of equations 
and plot of curves requires the use of a computer. 

We will proceed now to apply these theoretical considera- 
tions to the study of a practical problem, specifically, to 
the determination of the set of near-optimal parameters to be 
used in the operation of an engine-governor system. This 
problem was suggested by the reading of reference (4) and 
subsequent exchange of correspondence with one of its authors. 
In this paper, it is suggested that the determination of 


the parameter settings must be made by assuming that the 


dominant roots of the system occupy a. certain well-defined 
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position, specifically that the real and imaginary parts of 
“ке complex roots must have equal magnitude equal also to 
the magnitude of the third (real) root. 

The parameter space method was used to find new combina- 
tions of parameter values that would permit us to determine 
if the root locations, as assumed in the aforementioned 
paper, would correspond to an optimal solution, and to 
search for new solutions. In addition, the results permit 


study of related problems such as the problem of sensitivity. 
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III. PARAMETER PLANE STUDY OF AN ENGINE-GOVERNOR SYSTEM 


A. GENERAL CONSIDERATIONS 

Most gas engines exhibit a certain degree of dead time, 
that is, an interval of time between a change in fuel flow, 
and the corresponding change in driving torque. This time 
lag is made of a number of components of which the most 
important are: the delay in "sensing" the speed change that 
follows the load change; the delay within the governing 
mechanism during which the reaction of the sensing element 
is being translated into governor action; the time required 
for the governor action to correct the fuel injection pump 
setting; the delay that must take place before the corrected 
fuel change can be introduced into an engine cylinder; the 
time required to convert the corrected fuel change into 
engine torque. 

Dead time is found to be detrimental to stable operation 
and tends to increase ss ent errors. It also deter- 
mines the maximum rate at which the engine-governor combina- 
tion can be operated. The purpose of this study is to find 
the best isochronous governor parameters so that optimum 
transients, for the operation of gas engines of several 
types be obtained. 

The study of ideal governor settings in the presence of 
dead time will be made in the parameter space. This method 
leads to a computer process for analyzing the variable system 


Parameters. The process is restricted to a two parameter space. 
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B. ENGINE ASSUMPTIONS 
The driving torque incremental change AQ is related to 
the small variations of the governor actuator AZ by the 
equation: 
т E (9) 
-TS , | | 
where e is the dead time operator. 


The minimum amount of dead time T that can be expected 


is given by the following formulae: 


T = 20/N for a Diesel engine (10) 
T = 40/N for a two-cycle gas engine (G.E.)(11) 
T = 60/N for a four-cycle gas engine(G.E.)(12) 


in which N is the rotational speed in RPM. 

The assumptions made relative to the engine rotating 
parts are that the accelerating torque is equal to the dif- 
ference between the load torque and the driving torque as 


expressed by the following equation: 
ΠΟ. ΠΠ Ge) 


Combination of (9) and (13) gives: 


= 5 


TI u B 
39 3 AN = C¿ e AZ - AL (14) 


Change to a system of per unit variables gives the 
engine equation used in the analysis: 
ЕК5 


Ту Sn = e Z- L (15) 


C. GOVERNOR ASSUMPTIONS 
It is assumed that the isochronous pressure compensated 


type of governor is used. The differential equation 
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describing its dynamics is: 


2 


(57 + а15)2 = -u, (S + a,)n (16) 


where н. = (N,/2,)a,. 


D. ENGINE GOVERNOR COMBINATION 


Combination of equations (15) and (16) gives: 


DM πμ. 
EH x: 2, 229 5-193n e ml sea. )5. (17) 
1 m Hm T. 1 


E. USE OF THE PARAMETER SPACE METHOD 


We start with the engine-governor equation (Equation 17): 


Π., Шы N 
(S740, 87+ ae TS oy 2 I m = —S (Sta, )2 
M M M 





From equation (17), the C.E. can be extracted: 


И UA a 
53 + ҒЫ ЕЕ — аа. 2 A (18) 





Since we have three adjustable parameters in the charac- 
teristic equation, we can control the position of three of 
its roots. We assumed then that the system had three domi- 
nant roots similarly to what would happen if the delay in 
the same equation were equal to zero. For these three domi- 
nant roots we had two choices with respect to their intended 
location, that is, we could choose to have three real roots, 
or two complex conjugate roots and a real one. Because we 
Were interested in having a small rise time for the tran- 
sients of the system, we chose the second case, that is, we 


assumed that the system should exhibit a pair of complex 
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conjugate roots and a single real root close to the origin 
in the S plane. 
In accordance with this choice we made: 
ρωσ απο + gu 
We have then that: 


ТОО ен?) - (290) (19) 


5 2 


ο = (3040. - сЗ) + 1(3σ΄ω - w>) (20) 


: -TS 
The exponential term e becomes: 


SUIS Е a1 (-otjw) 


To =I OT 
e e 


(cos WT - j sin WT) (21) 


Substituting (19), (20) and (21) in (18) gives: 


Πο - o^) ЕСА = Ae: a, Co^ E 28 j20w) 


+ z^ e'?(cos wT-j sin wT) (-0+jw) 
M 
на 
+ c e "(cos wI-j sin wT) = 0 (22) 
M 


Both the real and the imaginary parts of this equation 


must be equal to zero: 


To 


u 
30W?-0°+u, (0?-u?)+ e (ш sin wT-0 coswT+ta_coswT)= 0 (23) 
M 


3 


2 3 H2 To ἃ ς 2 
30“w0-w*-20,00+ ш-е (0 sinuTtu cosuT-a,sinuT) = 0 (24) 


1 3 


Equations (23) and (24) form a system that can give solu- 
tions for any two parameters if the others are fixed. If we 
take a, and №, for parameters of the plane, and fix the value 


ik 


of о,» Че сап вес а complete representation of the loci of 
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the points (a, w), H, (w), С = constant or a, (0), μρίσ), 


constant] for which the C.E. has a solution of the form 


= 
I! 


-0 + ju. 

The curves determined as above-curves O and 0 аге loci 
of all parameter pairs (αι and μο) for which the solutions 
-o + 4W and -o-jw exist. 

If we denote by ô the dominant real root of the C.E., we 


can by the same methods get a plot of Š lines in the а,» 12 


parameter plane. 


We can repeat the process of finding o, w and 6 curves 


for as many values of E as we want. 


For a certain a, = constant, the intersection of selected 


5% апа wW, curves is the loci of values of the parameters 04 


and и, for which the selected solution 202 T jo, exists. 


In the three dimensional orthogonal space (a and 


JE 
а.) the loci of pairs aj. and и, is a line along which it is 


guaranteed that any readings of A» H3 and à, guarantee the 


3 : 
existence of the above form of solution. In a similar way, 
the par of lines of constant S» one for each value of о, 
forms a plane in the three dimensional space. 

The single or multiple intersections of the aforementioned 
line with the plane of constant 6$ - б о gives the set(s) of 
parameters 1» H3 and о; for which the system has a real 
root б and a pair of conjugate roots vU as ju ° However, 
since we are primarily interested in finding the near-optimal 


set of parameters for the governor, that is, the set of para- 


meters that offers a "better" transient, we will not be 
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looking for a set of the three parameter space parameters as 
indicated above. Instead, we will fix the location of the 
pair of complex roots, which, as explained above, will allow 
the determination of a line in the three dimensional para- 
meter space. We will then proceed to test the system with 
parameters corresponding to points of this line. Of course, 
the parameters corresponding to each of these points deter- 
mine a different location for the dominant real root. 

A computer program was written to plot ar 0, 0 and 6 
curves on a parameter plane with axis Ay and H. for a number 


of values of the parameter a, (computer program #1). This 


3 
Program was designed with the purpose of making it adaptable 
to studies of any type of engine-governor combination that 
may be described by equation (17). 

We performed the computer analysis first for a hypothe- 
tical system in which the transport lag was zero and later 
for each of the three practical cases considered--Diesel 
engine, two-cycle gas engine, and four-cycle gas engine-- 
each of them with the same C.E., but with different dead 
times as in equations (10), (11) and (12). The problem of 
obtaining parameter values that guarantee certain root loca- 
tions is thus solved. We still needed, however, any kind of 
way of testing the systems with each of the sets of para- 
meters obtained through the parameter space computer program. 
We, naturally, made use of a second computer program (compu- 


ter program #2) which was designed under the following heading. 
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IV. 


We want a 
ferent groups 


determine the 


TESTING OF THE ENGINE-GOVERNOR SYSTEM 


program to test each of the systems with dif- 
of parameter settings. From equation (17) we 


transfer function (T.F.): 





e S(S + a,) 
: Tu 1 
sae as + Le 8 + =." 
M M 
The denominator of the above T.F. is rearranged as 
follows: 
3 ж Yo -TS Ἡ 5  -τς 
DET 015 + m = S + T e 2 
М М 
2 ШЕ 15 
= S (S + Q.) + — e (S + a.) 
1 7 3 
M 
н, τ + а.) 
ы ο ο no 
M S (S + αι) 
Then, the T.F. becomes: 
- = S(S + αι) 
та = Ы 
x 2 н» аз ο + а.) 
ТЕСЕ - ) 
M S (S + a,) 
- ш-5 
" M р 
Т -TS ` l-GH (200 
u. e (S + 25) 
1+ = > - 
М S (S + o.) 
1 
ο тез πας: cuc (27) 
TS 
M 
ue ^ (Sa) 
H = (28) 
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Using block diagram description, we can partition G and 


H in the following form: 








And the complete block diagram of the feedback system becomes 
as represented in Figure 8. This configuration allows the 


use of C.S.M.P. 
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V. PROCESSING OF A PRACTICAL PROBLEM 


AND ANALYSIS OF RESULTS 


The design of ee requires the choice of a "stan- 
dard" type of response to load changes. Since a change in 
load unavoidably causes a Eo change in speed, the 
usual application has the following requirements: 

a) Recovery should be fast and with negligible differ- 

ence in steady state speed; 

b) The maximum deviation from steady state speed should 

be less than some specified percent; 

c) The transient recovery should not overshoot the 

steady state value. 
Trade-offs between these requirements may be necessary in 
some applications. 

For a system with no time delay, but with third order 
characteristics, the best solution for items a) and c) is 
a set of two complex and one real root, such that the real 
root is numerically equal to the real parts of the complex 
roots, and such that the imaginary and real parts of the 
complex roots are equal. For typical applications, a numeri- 
cal choice which gives a satisfactory en is a real 
root at -3 and complex roots at -3 + j3. (A study of re- 
sults based on other choices is made later.) Justification 
of this is given by Figures 9, 10, 11 and 12, which show the 
transient responses for a system with no time delay, with 


complex roots at -3 + j3 and with several locations of the 
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-— 


real root. Notice that Figure 11 is the one that corresponds 
to having the real root 6 equal to -3 and that Figures 9 and 
10 don't satisfy item a) while Figure 12 shows an overshoot 
which goes against the requirement posed in item c). 

We wanted to find the near-optimal parameter settings 
for the governor described in reference 4 for each of the 
types considered in the same paper. We assume in all cases 
that the steady state speed of the engine is 450 RPM. 

We used some practical data as it appears in the same 
reference: 

7.805788 


T = КЕ! 


20/450 for a Diesel engine 


407450 for a 2 cycle С.Е. 


I! 


60/450 for a 4 cycle G.E. 

We needed also to know what was the behavior of a system 
governed by equation (17), but with zero delay. This could 
be helpful in establishing the effect of the time delay and 
was used in determining parameter values that, entered as 
data in computer program #2, allowed the plotting of Figures 
JE 10, LI πα το. 

In each case the parameter plane program was run for 
values of the parameter о, eguat to 05, 0.6 0.7, 0.8; O 
EE 1 201.5, 1.4, 1.5, 1.55, 1.6 and 1.7. 

We decided in all cases that this system should have a 
pair of roots at -3 + j3 and so we concentrated on finding 
the location of the point of the parameter planes where the 


curves w = 3.0 and o = 3.0 intersected. Had the choice of 
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root locations been different, we could still have used the 
same parameter plane plots since they contain curves for 
ENS of c and o equal to 0.5, 1.0, 1.5, 2.0, 2.5, 3.0, 3.5, 
4.0, 4.5, 5.0, 5.5 and 6.0. The reason why we did not run 
the program for values of а, higher than 1.7 was that ме 


found that in each case for values of a, higher than a well- 


3 
determined value, the 9 = 3.0 and w = 3.0 curves on which we 
concentrated, would not cross in the parameter plane. This 


is justified by the following: 
For example, in the case of the four cycle gas engine 


with o = w = 3.0, we have from equation (23) that: 


54 + HW, x 0.17 x E sin 0.4-3 cos 0.4 + Q4cos 0.4)= O 
(254 Ы; (0.92 Q4 - 1.59)=-54 
2 920, - 1.59 


Assuming that the physical nature of the parameters re- 
quires that they have positive values has the effect that 
y, must be positive, which implies that 


1.59 
О το 





1.73 


The analysis of the parameter plane plots proceeded as 


const., we determined the intersection 


follows: for each о. 


of the oO = 3.0 and w = 3.0 curves and read the corresponding 


coordinate values of α1 


Plotting the point correspondent to these parameter 


and Ha: 


values in page 2 of the parameter plane plot, (page of curves 


of constant 6) allowed the determination of the location of 
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the real root for the considered value of a This value of 


3° 
ó was further checked by using computer program #3. 

Tables 1, 2, 3 and 4 list in each row the parameter 
values obtained by this process, and that correspondingly 
must be used in testing the dynamics of each of the systems 
so that they present a pair of complex roots at -3 + j3, 
and the value of the real root determined as described. 

. These tables were used to plot Figures 13, 14 and 15 
that represent the projections in each of the orthogonal 
planes that form the parameter space of the points of this 
Space that guarantee that each of the systems will have a 
ШЕШ ОГ roots at -3 + 33. Notice in these plots that the 
points correspondent to a certain system are disposed "in a 
line". Notice also that due to this form of disposition, 
if we were to work with a system described by equation (17), 
but with a delay different from those considered here, the 
determination of parameter values that would give the same 
root locations would be possible by simply using interpola- 
tion. 

Figures 16 through 55 are some of the parameter plane 
plots obtained by using the parameter plane program in the 
IBM 360/65 digital computer. These are only a part of the 
large number of parameter plane plots that we had to obtain 
to construct Tables 1, 2, 3 and 4. We chose those that 


either give parameter values close to those found to the 


"best" or the ones that give real "bad" parameter sets. 
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о ИЗО з 


145.8 -1.09 


158.0|-1.49 


9/5» 


IIb |-3.33 


242.0 |-3.882 


Table 1 - T = O 
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Table 2 


842.2 -0.60 


154.81|-1.15 


175.0 |М оз 


2354.0 |«2982 
238.6|-3.00 


2:5 3999-3 759 


306.0.) 028% 





T - 20/450 





145.8|-1.09 154.8|-1.15 


1150), 0 |-1.29 165.0|-1.29 


158.0|-1.49 


167.0 -1.68 


177.0 -1.95 


235.0|-2.88 


253.8 -3,45 


ATOR 0 37285 


306.0 |-ü4.88 





Table 1 - T = O Table 2 - T = 20/450 
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207.0) —0.90 


251.0 -1.25 


368.0 -2.15 





Table 3 - T = 40/450 Table 4 - T = 60/450 
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Each of the 3 systems was simulated including their res- 
pective time delays and the parameter values as in each row 
of tables 1-4. Figures 56 through 70 are plottings of the 
transients of each of the systems when ——— to 100% of 
load change. We obtained the transient responses for every 
set of parameters in those tables, but we are IEEE EE. 
as many as needed to illustrate the following analysis. 

Figure 71 is a plot of the transients of each of the 
systems for the case in which the dominant real root has a 
value of -3,0. We can see that all these transient res- 
ponses are very SOC ERES anas Thus the three roots considered 


in the cases of the systems in which the delay is non-zero, 


are really dominant over all the others, and the determina- 


tion of the near-optimal parameters for a system with trans- 
port lag consists of finding the parameter values that 
determine the Т of the dominant roots to be the 
ideal location of the roots for a system with the same 
equation, but with zero time delay. (This criteria is 


restricted to a certain portion of the S plane as will be 
ci later.) 

Examination of Figure 71 shows that for each of the 
systems, the settling time is very approximately the same; 
the differences in time delay determine only the amplitude 
of the transient, this being proportional to the amplitude 
of the delay, which makes sense since the greater the delay, 
the more time the governor takes to interpret the change 


and set the needed corrections. 
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A comparative analysis of Figures 56 through 70 shows 
that the transients of each system with parameter settings 
that correspond to having the real root closer to the origin, 
have larger settling times, while in the opposite case, the 
recoveries are faster, but an overshoot appears. 

The study of the method of determination of parameter 
sets is completed. We tried to make use of the large amount 
of data available by proceeding with the analysis of some 
other subjects. The first of these relates to the choice 
of location of the pair of «complex roots. 

For a linear system, we know that the same root pattern, 
that is, having the real root numerically equal to the real 
part of the complex roots and this real part equal in magni- 
tude to the imaginary part, should give the same type of 
transient, except of course for time scaling. In other 
words, the three groups of dominant roots, 

=з -4 -10 
э апа 
- 3558 =44+34 # 04310 
should all give identical transients for a linear system, 
— for the time scale. However, in the case of a system 
with time delay, this should not apply, and there should be 


some numerical limit on the pattern, depending on the value 


of T. In other words, if we tried to make this root pattern 
fall far to the left of the S plane, either the parameters 
would become negative-therefore physically unrealizable-or 
some other roots would become dominant, or would even fall 


in the right half of the S plane. The answer to this should 
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be looked at, because, given that this pattern were realiz- 
able more to the left of the S plane, eut mM tran- 
sient would become faster, and the magnitude of the maximum 
deviation could also possibly be smaller. 

We started by finding the sets of parameters that would 
give to the Diesel governor system root EEE ΤῊΝ 
at (-4 and -4 + j4), and € and -5 + j5). The determina- 
tion of parameter sets went as before, that is, in each case, 
we looked in page one of each of the plots of const. а, for 
the intersections of the O = á with ü = 4 and o = 5 with 
w= 5. For each of these intersections, we determined in 
page two of the plot the corresponding value of the real 
root 6, and were finally able to determine the values of 


the parameters о © апа Uo» that guarantee the desired 


T? 
root locations. The parameter sets determined are displayed 


in the following table: 













DOMINANT 
ROOTS 1 5 Ma 


3 
5, -5+j5 20.91 2225 643.0 


We used computer program #2 to determine the responses 












of the system while working with these sets of parameters, 
and we obtained the transients displayed in Figures 72 and 
73. Looking at Figure 72, it can be seen that, with the 
parameters set such that the system has dominant roots at 


-4 and -4 + j4, the speed transient is faster and of smaller 
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amplitude than any one of those previously found, which 
guarantees that the correspondent choice of parameters is 


"optimal". However, the response, displayed 


closer to the 
in Figure 73 of the system with dominant roots at -5 and 

5 + j5, is not good under item c) of the established cri- 
teria. Here the recovery is very fast and the amplitude 

of the transient is even smaller but an overshoot shows. 
This proves that our hypothesis about the existence of a 
limit in the process of bringing the root pattern more to 
the left of the S plane, was reasonable. 

Another topic of interest was the one related to toler- 
ances in the parameter values and/or sensitivity of the 
systems to parameter change. The possibility of and ease 
in performing such studies is one of the great advantages 
of the parameter space method, since the parameter plane 
plots display sets of curves that allow to relate changes 
in parameter values with changes in root location. Figures 
74 and 75 show a parameter plane plot where we marked the 
point that confers to the engine-governor system a set of 


roots at -A and -A + ДА. The -.- lines correspond to varia- 


tions of a certain percent in the determined values of the 


parameters а. апа н). The possible changes in root locations, 
due to changes in parameter values, is quite obvious. For 
example, the considered variation in the value of Qs while 


№, is kept constant, can cause variation of w between D and 
E and variations in the values of GO between B and C. 


Obviously, similar considerations apply to page 2 of each 
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plot and allow the study of changes in the value of 6. The 
analysis of the effects of changes in а, is a little more 
difficult, since it requires the use of plots correspondent 


to adjacent values of a Here we should plot the point 


3' 


correspondent to the nominal values of Q, and Uo in the 


pl 
adjacent plots. The point plotted in the adjacent plots 
shows immediately the changes in 0, O and Ó due to the varia- 
tion of 25 

Of course, if it were desirable to do so, the whole 
process could be reversed; that is, we could easily find 
out the new values of the parameters, if it were required 
to change the location of the roots in relation to the 
values set initially. 

Finally, notice that the parameter plane plots are a 
great aid in studies of sensitivity of the systems since 


the slope of the lines plotted gives direct measure of that 


eharaeteristäic, 
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VI. CONCLUSIONS 


The association of the parameter plane method, with 
appropriate computer analysis, converted a rather difficult 
problem into an almost elementary one. In effect, after 
passing the difficulties in writing and debugging the para- 
meter plane program, the analysis of the information avail- 
able is so easy, that an ill-advised observer may not notice 
all the trouble involved in the solution of the same problem, 
if it were attempted by other methods. 

The problem of finding the set of near optimal parameters 
with time delay consists of finding the parameter values 
that guarantee that the dominant roots of the C.E. of the 
system will be located according to a special pattern. The 
quality of the response improves with the displacement of 
the dominant roots to the left of the S plane but, this pro- 
cess has a limit depending on the value of the delay. The 
determination of this limit at which the quality of the res- 
ponse should be near-optimal has to be made through a trial 
and error process. 

The method described in this study is very easy to follov, 
and the information available allows for investigation in 
many areas of interest, such as sensitivity of the systems to 
variations in parameter values, analysis of permitted toler- 
ances in these values, optimization of performance under a 


defined cost function, and others. 
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Figure 1. Control systems with transport lag in the 
forward and backward paths. 
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Figure 2. Plot of simultaneous solution of equations 


E7) and (8). 


Figure 3. ος line in the À = const. plane. 
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Figure 7. M plane in A, B, C parameter space. 
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Figure 11. "Best" possible response of a hypothetical 
system with zero time delay. 
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with zero time delay. 
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Figure 59. Speed transient of Diesel. 


13 





ΓΕ 





Due 


mos 
--) 


o [y 


JUN 


-LIL 


Figure 60. 





б | 


i ; 
5 
! 


ea T= 20/4450 Аз-/7 


A1=/35 . MU2-2060 
X SCALE=0.50 UNITS/INCH 
Y SCALE=0.01 UNITS/INCH 


Speed transient of Diesel. 


94 





пас 


955 





-019 
Ba 


MN T= 40/450 Аз 05 

A172 25 MU>= /50.3 

X SCALE=0.50 UNITS/INUH 
"E Y SCALE=0.01 UNITS/INCH 


Figure 61. Speed transient of 2-cylinder G.E. 





pu 
--Δ- 


=i) 


-DZ 


Г 
% 
A 


QS 


ες: 


Γ 
` 


un 
c 


Figure 62. 





па 
ΠΗ 


A17 /46 MU? - 266.0 
X SCALE=0.50 UN TS/INCH 
Y SCALE=0.01 UNITS/INCH 


Speed transient of 2-cylinder G.E. 


96 





EN 


J 


e 
uk 


-032 


у^ 
=з 





‘O03 


Land 
— 3 


Т= 40/450 Аз- 134 


A1= 15.26 MÜ = 3014 
X SCALE=0.50 UNITS/INCH 
Y SCALE=0.01 UNITS/INCH 


Figure 63. Speed transient of 2-cylinder G.E. 


855 


97 





3 


09 


| | МИ | | 
Гу T= 40/450 A3=14 


A1= /6.1 MU 2 =320.4 
X SCALE=0.50 UNITS/INCH 
Y SCALE=0.01 UNITS/INCH 


sa s 
«Ὁ δ. 





-092 


m a 
ف اراس‎ 


~ 
-— 


una 
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